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ABSTRACT 

We explore the relationship between the models for progenitor systems of Type la supernovae and 
the properties of the supernova remnants that evolve after the explosion. Most models for Type la 
progenitors in the single degenerate scenario predict substantial outflows during the presupernova 
evolution. Expanding on previous work, we estimate the imprint of these outflows on the structure 
of the circumstellar medium at the time of the supernova explosion, and the effect that this modified 
circumstellar medium has on the evolution of the ensuing supernova remnant. We compare our 
simulations with the observational properties of known Type la supernova remnants in the Galaxy 
(Kepler, Tycho, SN 1006), the Large Magellanic Cloud (0509-67.5, 0519-69.0, N103B), and M31 
(SN 1885). We find that optically thick outflows from the white dwarf surface (sometimes known 
as 'accretion winds') with velocities above 200 km s -1 excavate large low-density cavities around the 
progenitors. Such large cavities are incompatible with the dynamics of the forward shock and the 
X-ray emission from the shocked ejecta in all the Type la remnants that we have examined. 
Subject headings: hydrodynamics — stars:binaries:close — supernova remnants — supernovae: general 
— X-rays:ISM 
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1. INTRODUCTION 

Despite decades of continuing effort, the progenitor 
systems of Type la super novae (SNe) have not been con- 
fidently identified yet fsee lHillebrandt &: Niemeverll2000l . 
and references therein). This has become one of the most 
pressing problems in stellar evolution, with profound im- 
plications ijo^_rosjnology aJidthe chemical evolution of 
galaxies (Kobavas hi et al.l [T998L It has been known for 
some time that Type la SNe are the result of the ther- 
monuclear explosion of a C+O white dwarf (WD) that 
is destabilized by accreti on of matter from a c ompanion 
in a close binary system ([Hovle &: Fowlerlll960f ). but the 
details of this process are still the subject of heated de- 
bate. In particular, the nature of the companion and the 
physical processes involved in the evolution of the binary 
prior to explosion remain obscure. 

Most of the discussions on Type la progenitor systems 
are held in the framework of two basic scenarios. In the 
so-called single degenerate (SD) scenario, the WD com- 
panion is a normal or slightly evolved star. In the double 
degenerate (DD) scenario, it is another WD. In the SD 
scenario, the WD might explode either close to the Chan- 
drasekhar limit (SD-Ch) or some time before reaching it 
(sub-Chandrasekhar explosions, SD-subCh). Several ob- 
servational and theoretical arguments point at SD-Ch 
explosions as the most promising route to Type la SNe, 
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but there are still many uncert ainties associated with this 
parad igm (for discussions, see iBranch et al.lll995t iToutl 
2005). The main problem is the necessity to bring the 
mass of the C+O WD close enough to the Chandrasekhar 
limit to trigger the explosion. At the onset of the fi- 
nal mass transfer episode in SD sys tems, the WD mass 
is thought to lie between 0.7 M© (Lan ger et al.l l200"(| 
lHan fc Podsiadlowskil[200l and 1.2 M g (|Nomotolll982( ). 
This means that, in order to achieve the ~ 1.38 M nec- 
essary for a thermonuclear runaway, the WD needs to 
accrete between 0.2 and 0.7 M G of material from its com- 
panion and then burn it somehow to C and O avoiding 
the highly unstable degenerate conditions at the WD sur- 
face. Our present understanding of Type la SNe in the 
SD-Ch scenario is based on theoretical models that at- 
tempt to explain the details of these accretion and burn- 
ing processes. 

One key aspect of such models is the presence or ab- 
sence of outflows from the binary system during its pre- 
SN evolution. If these outflows are indeed present, they 
should leave some kind of imprint on the circumstellar 
medium (CSM) of the progenitor, but the evidence for 
this modified GSM in Type la SN observations remains 
inconclusive. Prompt emission at radio or X-ray wave- 
lengths from the interaction of SN ejecta with surround- 
ing material, which has been found in several core col- 
lapse SNe, appears to be ab sent in Type la events. In the 
radio. IPanagia et ail (|2006f ) report no VLA detections in 
27 nearby objects. In the X-rays, no Type la SNe have 
been detected with RO SAT (upper limits on one object, 
iSchlegel &: Petrell993f ) or Chandra (upper l imits on three 
object s, Hughes et al. in preparation). Ilmmler et al.l 
(2006) report Swift observations of eight Type la SNe, 
with only one marginal (~ 3a) detection, SN 2005ke, a 
peculiar object that was not seen by Chandra. It is im- 
portant to keep in mind that these studies only probe 
the CSM at distances between 10 15 cm (for observations 
weeks or months after the explosion) and 10 17 cm (for the 
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latest radio observations in iPanagia et al.l 120061 . several 
years after the explosion). Modification of the CSM on 
larger scales will become apparent in the aftermath of the 
ex plosion, during the supe rnova remnant (SNR) phase. 
In lBadenes fc Bravol (|2001l ) (henceforth BB01), we con- 
ducted a preliminary search for the imprint of presuper- 
nova outflows on the dynamics of young Type la SNRs. 
In the present work, we expand and revise the calcula- 
tions presented in BB01, and we perform a more detailed 
and robust comparison to the observations of SNRs. Our 
aim is to use these comparisons as a means to constrain 
the fundamental properties of the presupernova outflows 
and gain insights on the nature of the progenitor systems. 

This paper is organized as follows. First, we examine 
the outflow mechanisms from Type la progenitors in the 
SD-Ch channel, and we build simple models to mimic 
them (§0). From these outflow models, we simulate the 
structure of the CSM around the progenitor system at 
the time of the SN explosion (§ [3]) and the subsequent 
evolution of the SNR (§[3]). In §[5j wc compare our SNR 
models with the known properties of the forward shock 
dynamics and X-ray emission in a sample of known Type 
la SNRs. Finally, we discuss our results in § [5] and we 
summarize them in § [7J 

2. TYPE IA PROGENITORS IN THE SINGLE 
DEGENERATE CHANNEL: PRESUPERNOVA 
OUTFLOWS 

2.1. Pre-SN Evolution Scenarios 

There are two channels by which a SD binary sys- 
tem can lead to a Type la SN. The systems in the first 
channel consist of a WD and a main sequence (MS) 
or slightly evolved companion star (WD +MS c hannel , 
iLi fe van den Heuvell lHachisu et al.lll999bl) . Mass 
transfer occurs via Roche lobe overflow, which requires 
that the two components be very close, with an initial 
period of the order of days. A different class of Type la 
progenitors systems is composed of a WD and a red giant 
(RG) star with a wider separation and a longer initial pe- 
riod, of the order of a hundred days (WD+RG c hannel , 
ILi fc van den Heuvell Il997t lHachisu et ail Il999af) . The 
mechanism for mass transfer in these more detached sys- 
tems is not well understood, but Roche lobe overflow 
might still play a significant role. 

The details of the evolution of Type la progenitor sys- 
tems in these two channels are very complex. This is 
mainly due to the fact that the process of mass trans- 
fer affects both the orbital parameters of the system and 
the evolution of the components, which in turn modify 
the conditions for the mass transfer. In the most widely 
accepted picture for SD Type la progenitors, the H- or 
He-rich material from the donor star forms an accretion 
disk and then accumulates in a shell on top of the WD. 
Given the appropriate conditions, nuclear burning in this 
shell converts the accreted material to C and O, until 
the mass of the WD reaches ~ 1.38 M . At this point, 
a runaway thermonuclear flame ignites in the central re- 
gions of the WD, and a Type la SN explosion ensues. 
It is important to stress that a SN explosion is only one 
of the many possible outcomes of the evolution of close 
binary systems that contain a WD. The key parameter 
that determines whether a specific binary system leads 
to a Type la SN or not is the accretion rate of the WD, 
M. 



For accretion rates below a critical limit M sta Mei nu- 
clear burning on the surface of the degenerate WD is 
unstable (|Nomotolll982D 6 . The value of M s tabie dep ends 
on the mass of the WD (|Nomoto fc Kon do 19 91J), on 
the m etallicity of the accreted material dPiersanti et al 



2000) , and on the rota tion of the WD (jUenishi et al 



2003tfYoon et al.ll2004al lbT), but it is generally assumed 
to be of the order of 10 -7 M yr _1 . If M is only slightly 
lower than M stable, nuclear burning proceeds cyclically 
in weak flashes, and most of th e accumulated materia l 
might be retained by the WD (|Kato fe Hachisul 12004?) . 
Recurrent novae like RS Oph or U Sco are thought to 
be the observational counte rparts of white dwarf bina- 
ries accreting in this regime (lHachisu fe Katoll200lh . For 
even lower values of M, however, the flashes become 
stronger, effecti vely turning into nova explosions that 
erode the WD (jYaron et alj 120051 ). If the accretion en- 
ters this regime, the WD mass effectively decreases with 
time and t he system cannot produce a Type la SN (see 
Fig. 5 in iKahabka fe van den Heuvell 119971 and refer- 
ences therein). 

At higher accretion rates, the energy input from 
the shell nuclear burning can lead to the expansion 
of the outer layers of the W D to red giant dimen- 
sions (|Nomoto fe Kondol 11991 . This occurs when 
M approaches the nuclear Eddington-accretion rate 
M E dd,nuc = L Edd /e nuc , where L E dd = AttcGM/k is 
the Eddington luminosity of the WD, and £ nuc ~ 7 x 
10 18 erg g _1 is the energy gained by fusing H into He. For 



a WD of one solar mass, MEdd.nuc — 6 x 1O _7 M yr" 
Upon expansion of the WD, the binary system may en- 
ter a common envelope phase, which might lead to a 
spiral-in (due to d ynamical friction) and the merging of 
both components (|Iben fc Livid 119931 : llvanova fc Taaml 
2004). For a long time, it was thought that this would 
prevent man y SD systems from b ecoming Type la pro- 
genitors, but lHachisu et all (|1996f) (henceforth, HKN96) 
proposed a mechanism to avoid the common envelope 
phase. These authors claimed that, for high values of M, 
the luminosity from the shell nuclear burning can drive 
an optically thick 'accretion wind' outflow from the WD 
surface with an outflow mass loss rate M f such that 
the effective value of M remains below M-Edd,nuc- This 
accretion wind results in the loss of material from the 
binary system. 

2.2. Accretion Wind Outflows from Type la Progenitor 

Systems 

In recent times, most of the models for SD Type la 
SN progenitors have relied to a greater or lesser extent 
on the accretion wind mechanism of HKN96 as a means 
to limit the effective value of M . This is despite the 
fact that several fundamental aspects of the accretion 
wind mechanism are not well understood. In particular, 
some kind of self-regulation is required so that M falls 
below MEdd,nuc (to avoid the expansion of the WD), 
but at the same time remains very close to this criti- 
cal value (in order for the WD to gain enough mass to 
explode as a Type la SN). Furthermore, the accretion 



6 In a recent paper, Starrfield et al. (2004) proposed a stable 
accretion regime for low values of M, but see Nomoto et al. (2006) 
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wind is supposed to remain active up to mass trans- 
fer rates of the order of the Eddington accretion rate 
M E dd,acc = L E dd/£acc, where e acc = GM/R is the spe- 
cific potential energy gained by accreting matter onto a 
WD with radius R. In this situation, M t ~ MEdd,acc, 
which implies, for wind velocities similar to the WD es- 
cape velocity (see § I2.2.1j) , a ratio of wind-to-photon mo- 
mentum of T] = l/cy/2Rn/{GM) ~ 10 2 . It is known 
that winds driven b y photon scattering i n lines cannot 
work in this regime (jLucv fc Abbottl fl993). In principle, 
continuum-driven winds might reac h such high efficien- 
cies, as claimed by HK N96 (see also lOwocki et al.ll2004t 
ISmith fc Owockil 12006). but this has never been proved 
with detailed simulations. Due to these and other con- 
cerns, the role that accretion winds play in the presu- 
pernova evolution of Type la progenitors remains con- 
troversial. Some authors claim that this mechanism can- 
not work in a realistic scenario, an d that H-accreting 
WDs cannot e xplode as Ty pe la SNe (jCassisi et al.lll998t 
iPiersanti et~aT1ll999L 12000 ). 

2.2.1. Theoretical studies 

Following the work of HKN96, lHachisu et all (|1999bf ) 
explored several specific cases within the WD+MS chan- 
nel, making use of analytical formulae to follow the evo- 
lution of the binary system pa rameters and mass tran sfer 
rates. In a companion paper, lHachisu et alj (|1999al ) ap- 
plied the same techniques to systems in the WD+RG 
channel. In these papers, substantial outflows from the 
Type la progenitors were found in all cases. The typi- 
cal timescales for the outflows are several times 10 5 yr, 
with mass loss rates M a f of the order of 10~ 6 M Q yr -1 . 
In some binary systems, the WD explodes as a Type la 
SN while the optically thick accretion wind is still ac- 
tive. In other cases, the outflow stops some time before 
the explosion, leading to a more or less extended mass- 
conservative phase. 

Other studies have followed the evolution of binary 
systems with accretion winds, relying on stellar evolu- 
tion codes ins t ead of the analytical formulae used by 
lHachisu et al.l (1999b. a). This kind of calculation is 
necessary in order to estimate the mass-transfer rates 
in the system correctly. In general, stellar evolution 
codes predict lower values of M f and hence longer evo- 
lutionary timescales than those obtained using analyt- 
ical formulae. The first s uch study was published by 
iLi fe van den Heuvell (|1997l ). but a more detailed and up- 
dated exploration of the paramet er space (with ~ 2300 
binary systems) can be found in lHan fc Pod siadlowski 
(2004). These works have in common that they found 
substantial outflow s from the Type la progenitors in all 
cases. In contrast, iLanger et al.l (j2000l ) only found out- 
flows in about one third of the progenitor systems of solar 
metallicity that they simulated. This is due to more op- 
timistic assumptions concerning the mass accumulation 
efficiency in accreting WDs, which might be justified by 
model s that include the effect of rotatio n (lYoon et al. 
l2004a|) . It is also important to mention that ILanger et al" 
stopped their simulations for systems which required out- 
flows with M a f > 3M E dd,nuc- 

In the theoretical framework laid down by HKN96, 
accretion winds are formed in the WD envelope with 
photospheric velocities of ~ 1000 km s -1 (Hachisu et al. 



1999b. a). These high outflow velocities are necessary be- 
cause the escape velocity from the surface of the WD 
is also very high (roughly of the same order of mag- 
nitude, see Figure 4 and accompanying discussion in 
iKato fc Hachisu! [19941 ). The orbital velocities in Type 
la progenitor syste ms are much lower than this, of the 
order of 200kms _1 (jLanger et al.ll2000D . 

2.2.2. Observations 

It is useful to compare the theoretical accretion wind 
models described in the previous section to the observa- 
tions of binary systems that have been proposed as Type 
la progenitor candidates. The list of candidate objects 
is very heterogeneous, and includes tr ansient supersoft 
X-ray sources like RX J0513.9-6951 ([Hachisu fc Katol 
2003b) , eclipsing bi naries with intermittent X-ray emis- 
sion like V Sagittae (Hachisu fc Kato 2003a) , and recur- 
rent novae like RS Oph (|Hachisu fc Katol 120011 ) and U 
Sco (jHachisu et al.ll2000D . 

Among these objects, RXJ0513. 9-6951 and V Sagit- 
tae are of particular interest, because they have non- 
explosive outflows that c an be modeled using t he ac- 
cretion wind mechanism (|Hachisu fc Kat o 2003b a|). In 
the case of RXJ0513. 9-6951, the outflows have a ve- 
locity of ~ 3 800 km s" 1 , with strong hints of a bipo - 
lar structure (|Pakull et al.lll993t iHutchings et al.ll2002j) . 
The outflows from V Sagittae are more complex, but a 
high velocity c omponent at ~ 1500 km s" 1 also appears 
to be present (|Wood fc Loc klev 2000). These high out- 
flow velocities support the theoretical values discussed 
in the previous Section. It is interesting to note that 
both systems are variable sources at optical and X-ray 
wavelength s, with periods of the o rder of a few hun- 
dred days. lHachisu fc Kato (2003b) attribute the vari- 
ability of RXJ0513. 9-6951 to the impact of the strong 
(~ 2 x 10 -6 M Q yr -1 ) accretion winds from the WD at- 
mosphere on the surface of the donor star, whose outer 
envelope is stripped off, temporarily shutting down the 
accretion process and interrupting the accretion wind 
itself. The same authors propose a somewhat similar 
model for V Sagittae, but it is unclear whether this kind 
of variability should be ubiquitous among binary systems 
with accretion winds. 

2.3. Modeling the Outflows 

Based on the theoretical studies described in § 12.2.11 
we have produced several models for the outflows from 
Type la SN progenitor systems. The temporal evolution 
of the outflow mass loss rate M j is shown in Figure [TJ 
and the fundamental properties of the models are listed 
in Table [TJ Model HI mimics the mass loss of the exam - 
ple system discussed in § 3.5 of lHachisu et al.l (|1999bf ) 
(see their Fig . 7). Model LV1 is built aft er the example 
in Figure 1 of lLi fc van den Heuvell (| 19971 ). Models HP1, 
HP2, and HP3 are approximations to the outflows from 
the three represe ntative Type la progen i tor sy stems dis- 
cussed in § 3 of lHan fc Podsiadlowskil (|2004h (Fig. la 
for model HP1, Fig. lc for model HP2, and Fig. le 
for model HP3). Finally, models LI and L2 correspond 
to the binary syste ms numbers 2 and 31 in Table 2 of 
ILanger et ail (l2000f) (see Fig. 7 for model LI and Fig. 
34 in iDeutschm ann 1991, for model L2). 

These outflow models are chosen as representative ex- 
amples of the typical evolution of Type la progenitors of 
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solar metallicity. With them, we expand and improve on 
the work presented in BB01, aiming at a more complete 
exploration of the parameter space for Type la progeni- 
tor systems. Our outflow models are evidently very sim- 
plified, but we believe that they capture the essence of 
the relationship established between the current models 
for Type la progenitor systems and the CSM (a more 
detailed discussion of some of the issues concerning the 
outflow models is deferred to § 16. ip . We note that, al- 
though we have sampled binary system calculations with 
different initial conditions and from different authors, the 
time scales and mass loss rates of all the outflows are 
very similar, and should shape the CSM and influence 
the evolution of the SNRs in a similar way. 

3. THE SHAPING OF THE CIRCUMSTELLAR MEDIUM 
3.1. Method and Parameters 

We have performed hydrodynamic simulations to cal- 
culate the impact that the outflow models presented in 
§ 12.31 would have on the structure of the CSM around 
the Type la progenitors at the time of the SN explo- 
sion. For this purpose, we have used VH-1, a numeri- 
cal hydrodynamics code developed at the University of 
Virginia by John Hawley, John Blondin and collabora- 
to rs. VH-1 is based on t he pie cewise parabolic method 
of IColella fc Woodward! (|1984h . and the code has been 
extensively tested and validated in multiple astrophysi- 
cal scenarios. For the timescales and velocities involved 
in the outflows from Type la SN progenitors, radia- 
tive cooling is dynamically important. We have in- 
cl uded this process thro u gh th e standard cooling curves 
of ISutherland fc Dopital (119931) for temperatures ab ove 
10 4 K, and the values of iDalgarno fe McCravl (|1972f ) for 
10 2 < T < 10 4 K (fo r a discussion of the cooling curves, 
see lRaga et al.lfl997h . The simulations are performed in 
one dimension with spherical symmetry, and they do not 
include thermal conduction. 

Once a particular outflow model is selected, the only 
parameters that need to be adjusted in our calculations 
are the outflow velocity u f, and the density and temper- 
ature of the interstellar medium (ISM), pism and Tism- 
Following the discussion in § 12.2.11 and § 12.2.21 we have 
chosen a reference value of 1000 km s -1 for u a f. For the 
ISM parameters pism and Tism we have chosen refer- 
ence values typical of the warm atomic ph ase of the ISM : 
p ISM = 10- 24 gcm- 3 and T 7SM = 10 4 K (lFerrierell200fl ). 
After a brief theoretical discussion (§ 13. 2[) . we will present 
the CSM structures generated by the outflow models us- 
ing these fiducial values (§ 13. 3| . and then we will examine 
the consequences of varying them (§ 13. 4j) . 

3.2. Theory 

A complete theoretical overview of the interaction be- 
tw een stellar winds and th e surrounding medium is given 
bv lKoo fc McKee] l| 1992a.b'). who built on the pioneerin g 
work of ICastor et alj (| 19751) and I Weaver et all (|1977l ). 
In the framework set by this theory, we can assimilate 
the accretion wind outflows to stellar winds, and expect 
them to form cavities or 'accretion wind-blown bubbles' 
around the Type la progenitors, similar to the more 
familiar stellar wind-blown bubbles found around the 
mass ive progenitors of core-collapse SNe (see iDwarkadasI 
120051 and references therein). In general, these cavities 



will have the usual structure of two shocks and a con- 
tact discontinuity (CD) separating four regions from the 
inside out: freely expanding outflow, shocked outflow, 
shocked ISM, and unperturbed ISM. The detailed con- 
figuration and radial extent of the cavity will depend 
on the temporal evolution of the outflow mechanical lu- 
minosity L Q f — {\/2)M fU^p the duration of the mass 
loss episode, and the pressure exerted by the ISM (for a 
review of all the possibilities, see § 7 in IKoo fc McKeel 
1992b). One parameter that is fundamental for the evo- 
lution of the CSM structure is the critical outflow velocity 
u cr , given by 



10 4 



M, 



°f u of Pism 
Ph 



l/n 



(1) 



where pu — 2.34 x 10~ 24 g is the mean mass per H 
atom in a solar abundance gas and all magnitudes are 
in c.g.s. units (expr ession adapted from eq. 2.5 in 
IKoo fc McKeel fl992a). Outflows with u a f > u cr are 
'fast', meaning that radiative losses do not affect the 
shocked outflow, and the expansion of the outer shock 
is driven by the thermal energy of the shocked mate- 
rial inside the cavity (energy-driven bubbles). Outflows 
with u a f < u cr are 'slow', and radiative losses affect the 
shocked outflow, either at the reverse shock or at the CD. 
These cavities are driven by the ram pressure of the out- 
flow itself (momentum-driven bubbles). For the values 

of M a f and Pism typical of accretion winds expanding 
into the warm ISM, the outflows are in the fast regime 
at velocities above ~ 100 km s -1 (see Figure [2]). 

3.3. Results: Reference Values 

The structure of the CSM at the time of the SN ex- 
plosion for each outflow model at the reference values of 
u of, Pism, and Tism is shown in Figure [3l The funda- 
mental properties of the CSM profiles are summarized 
in Table [2] As expected for outflows in the fast regime 
with timescales of ~ 10 6 yr, all the accretion winds ex- 
cavate very large low-density cavities around the Type 
la progenitor systems, with radii R c between 17 and 36 
pc. The temporal evolution of these cavities can be di- 
vided into two distinct phases. During the first phase, 
the outer shock is radiative, and the bubble expands at 
supersonic velocities. Thin shells of radiatively cooled 
material develop behind the outer shock, and then be- 
come thicker as the temperature of the material drops to 
10 2 K and cooling ceases. In the terminology introduced 
bv IKoo fc McKeel (fl992bl ). this is an adiabatic bubble 
with a radiative outer shock (ABROS), but the only ac- 
cretion wind-blown cavity that is still in this stage at ts n 
is that generated model HP3. Eventually, the pressure 
inside the cavity drops below the ISM pressure (which 
is Pism = PismTism^Ip = 1-33 x 10~ 12 dyncm~ 2 for 
the reference ISM parameters), and the cavity becomes 
pressure-confined. During this stage, the outer shock 
becomes a standard sound wave that rapidly dissipates, 
and the bubble expands (or contracts) at subsonic speeds 
(adiabatic pressure co nfined bubble, APCB, see § 4.4 in 
IKoo fc McKeel I1992H ). Models HI, LV1, LI, L2, HP1, 
and HP2 are in this stage at tsN- The reverse shock 
can only be seen in the two models that have an active 
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outflow at tsN (HP3 and LI). In these cases, the progen- 
itor system is surrounded by a small region of unshockcd 
outflow with mass M uo f that extends to a radius R uo f ■ 

3.4. Results: Modified Values 

Among the seven outflow models presented in Table 
[TJ we have chosen models HP3 and L2 as representative 
examples of outflows producing 'small' and 'large' cavi- 
ties. Based on these models, we have explored the range 
of CSM structures that might be found around Type la 
progenitors by modifying the properties of the ISM and 
the values of u t. 

Variation of the ISM properties. — The reference values 
of pism and Tism are at the high end of the expected 
range s for the warm atomic phase of the ISM (|Ferrierel 
2001). This choice maximizes the value of pism, so 
that the pressure containment of the wind-blown cavi- 
ties is maximum, and the size of the wind-blown bubbles 
shown in Figure [3] is a lower limit. To evaluate the ef- 
fect of decreasing pjs m to the lowest reasonable value for 
the warm ISM, we have evolved the outflow models HP3 
and L2 into an ISM with pism = 5 x 10~ 25 g cm~ 3 and 
Tism — 5 x 10 3 K. This translates into an ISM pressure 
of pism — 3.33 x 10~ 13 dyncm -2 , roughly an order of 
magnitude below the reference value. The results can be 
seen in Tableland Figure [4] labeled as models HP31owp 
and L21owp. These bubbles are ~ 20% larger than the 
reference cases, but there are no significant differences in 
the structure of the cavities, because the outflows stay in 
the fast regime. Even at the lower value of pism, model 
L21owp is already making the transition to the APCB 
stage at tsN- 

Variation ofu a j. — We have calculated the CSM struc- 
tures obtained by models HP3 and L2 with outflow ve- 
locities of lOkms -1 (a reasonable lower limit, being the 
typical value for the sound speed in the ISM), 100 km s -1 , 
and 200 km s -1 . The results of these calculations are 
presented in Table [2] and Figure [5] (models Hp3ule6, 
Hp3ule7, Hp3u2e7, L2ule6, L2ule7, and L2u2e7). The 
lowest values of u f in this sequence (10 and 100 km s -1 ) 
lead to outflows in the slow regime and momentum- 
driven cavities with structures that are very different 
from the ones we ha ye seen up to now. T he types of 
cavities described by iKoo &: McKed (|1992al) are nicely 
laid out by the sequence of HP3 models with increas- 
ing values of u f. Model HP3ule6 is a radiative bubble 
(RB), with a wind shock that is still radiative at tsN', 
model HP3ule7 is a partially radiative bubble (PRB), 
with an adiabatic wind shock, but active cooling of the 
shocked outflow at the CD; model HP3u2e7 is already 
in the fast regime, and has become pressure confined at 
tsN (APCB). The sequence of L2 models reveals an in- 
teresting effect associated with accretion winds that have 
a mass conservative phase prior to the SN explosion. If 
these outflows are in the fast regime, the energy-driven 
cavity survives for a long time, although it may become 
pressure confined and start to shrink (c.f. model L2u2c7, 
which is an APCB at is at). If the outflow is in the slow 
regime, on the other hand, the momentum-driven cavi- 
ties collapse when the accretion wind ceases. In the case 
of model L2ule6, which is a RB during the active out- 
flow phase, a relic shell of cooled material is left behind, 



but the PRB of model L2ule7 (which looked very similar 
to HP3ule7 in the active phase) has disappeared almost 
completely at tsN- 

We conclude this section with a brief reference to our 
previous results from BB01. The hydrodynamic simu- 
lations in that paper were very simplified, and did not 
include radiative losses. The cavity sizes for the outflows 
in the fast regime (models A and C in BB01) are quali- 
tatively correct, but the results for models with outflow 
velocities of 20kms _1 (models B and D), which are in 
the slow regime, are incorrect (the simulations were later 
repea ted with the inclusion of radiative losses in lBadenesI 
l200l . The main conclusion of BB01 - that fast accre- 
tion wind outflows lead to large cavities around Type la 
progenitors - still holds, but the more detailed results 
that we present here supersede the calculations in our 
previous work. 

4. SNR DYNAMICS IN THE MODIFIED CIRCUMSTELLAR 
MEDIUM 

The dynamics of SNRs expanding into low-density 
wind-blown cavities has been studied by several a uthors 
(|Tenorio-Tagle et al.lll990l 119911 : lDw"arkadasll2005t ). The 
reader is encouraged to consult these works for discus- 
sions - we will not go into the details here. Essentially, 
the SN ejecta expand almost freely until the inner edge of 
the cavity R c is encountered. What happens after that 
depends on the momentum of the ejecta and the mass 
that is co ntained in the she ll that surrounds the cav- 
ity, M s h- iDwarkadasI (|2005l ) described this interaction 
in terms of the A parameter, defined as the quotient be- 
tween the shell mass and the ejecta mass, A = M~ s h/M e j. 
Among the CSM structures presented in § O those gen- 
erated by outflows in the fast regime always correspond 
to the more extreme A >> 1 scenario (c.f. values of M s h 
listed in Table [2]). In this case, the forward shock (FS) 
becomes radiative when it hits the shell, and the Sedov 
stage (i.e., R fs oc t 0A ) is absen t from the SNR evolution 
(see § 7.3 in IDwarkadasI [20051) . This strong modifica- 
tion of the dynamics will have two important observa- 
tional consequences. First, the fundamental properties 
of the FS (radius, velocity, and expansion parameter) at 
a given SNR age will be very different from the case of 
an interaction with a uniform ISM. Second, since most of 
the ejecta will expand to very low densities before being 
overrun by the reverse shock, the X-ray emission from 
the shocked ejecta in the SNR will also be affected, lead- 
ing to a much lower ionization state at a given age. The 
CSM structures generated by outflows in the slow regime 
are more diverse, and their imprint on the FS properties 
and the X-ray emission from the shocked ejecta will not 
be so dramatic. 

In order to produce SNR models that can be compared 
with the observations, we have used VH-1 to simulate the 
interaction of t he Ty pe la explosion model PDDe from 
iBadenes et al.l (|2003l ) with several of the CSM profiles 
presented in § O The choice of a particular Type la SN 
explosion model does not have a significant influence on 
the radius and velocity of the FS, as long as the kinetic 
energy produ ced in the explosion is ~ 10 51 erg (see Figs. 
3a and 3c in IBadenes e t al. 2003]). The X-ray emission 
from the shocked ejecta, on the other hand, is profoundly 
affected by the differences in density and chemical com- 
position profiles of Type la SN explosion models. We 
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have chosen model PDDe because it maximizes the ion- 
ization state of the elements in the shocked ejecta for a 
given value of pism (compare panel 6d with the other 
panels of Figure 6 in lBadenes et al.| [2003). We have fol- 
lowed the noncquilibrium ionization (NEI) processes in 
the shocked ej e cta w ith the ionization code described in 
iBadenes et al.l ([20031 ). with one important modification. 
Some extreme cases of CSM interaction that we will spec- 
ify in § !5.3l lead to catastrophic cooling in the outer ejecta 
layers. To model this process, we have included ioniza- 
tion and radiative losses in the code, using the latest 
atomic data pub l ished by the CHI ANTI collaboration 
(|Dere et al.lll997t lLandi et al.ll2006l ). To avoid recalcu- 
lating the hydrodynamic evolution in each case, we have 
assumed isobaric cooling in the affected layers (i.e., the 
layers cool but remain in hydrostatic equilibrium with 
their surroundings), which is an excellent approximation 
for the inter-shock region in one-dimensional hydrody- 
namics. 

5. COMPARISON WITH OBSERVATIONS 

5.1. Observations: Forward Shock Dynamics and X-ray 
Emission from the Shocked Ejecta 

We have compiled the relevant observational data for 
young SNRs that have a firm Type la identification in Ta- 
bles[S]and|U We have restricted ourselves to objects that 
are close enough to study the dynamics of the FS, and 
(if possible) the X-ray emission from the shocked ejecta, 
and whose age is either known (i.e., historical SNRs) or 
can be estimated with a high degree of confidence. This 
yields seven objects: SN 1885 in M31; Tycho, Kepler and 
SN 1006 in our Galaxy; and 0509-67.5, 0519-69.0, and 
N103B in the Large Magellanic Cloud (LMC). Among 
these objects, SN 1885 was identified as a Type la SN 
from the optical spectrum recorded in the nineteenth cen- 
tury (|de Vaucouleurs fe Corwinlll985f). although it was 
probably a sublu minous event ([Chevalier fe Plai^ Il988t 
iFesen et al.ll2006l ). Tycho and SN 1006 have been tradi- 
tionally regarded as prototypical Type la SNRs based 
on their X-ray spectra and other evidence, and their 
origin has been confirmed by detailed modeling of the 
ejecta emission (for Tycho, see Ba denes et af] 120061 for 
SN 1006, Badenes et al., in preparation). The only 
object in our sample whose origin is somewhat contro- 
versial is the Kepler SNR (for a disc ussion and refer- 
ences, see ICassam-Chenai et al.l I2004T ) . but the promi- 
nent Fe emission and the virtual absence of O in the 
shocked eje cta revealed by the r ecent deep Chandra ob- 
servations (|Revnolds et al.l [2 006[ ) are strongly indicative 
of a Type la origin. The LMC SNRs 0509-67.5, 0519- 
69.0, and N103B w ere first identified as Type la by 
iHughes et al.1 (|1995l ). The results of later works appear 
to ag ree with the Type la hyp othesis in the case of 0509- 
67.5 (jWarren fc Hughes! l2004ft and N103B ([Lewis et all 
I2003h . but there are no detailed models in the literature 
that can confirm this. An interesting development for 
these LMC SNRs has been the recent detection of light 
echoes from the explosion, which has greatly reduced the 
uncertainty in their ages and might eventually confirm 
their Type la origin ([Rest et al.ll2005D . 

The radius and velocity of the FS in a SNR provide 
a first assessment of its dynamical state. The FS radii 
Rfs are determined from the angular radii ckfs and the 
distances to the SNRs D (Rfs = ctFsD). All the angu- 



lar radii listed in Table [3] are from Chandra observations 
except that of SN 1885, whic h is from the abs orption 
features seen in HST images (|Fesen et al.l [2006). The 
X-ray radii are very precise: even in the LMC SNRs, 
the Chandra PSF will only introduce errors below ~ 5%. 
The value o f D is know n with some accu racy for the 
LMC SNRs (|Alvesll2004D and for SN 1885 ([Fesen et all 
2006), but it can be very uncertain in the Galactic SNRs. 
The FS velocities are obtained from the study of the 
narrow and broad c omponents of the Ha emiss ion in 
non-radiative shocks (|Chevalier fc Ravmondl [l978). The 
derivation of ufs from the width of the broad Ha com- 
ponent is not straightforward, and it requires some mod- 
eling of the plasma physics at the shock. It is important 
to note that the FS velocities are obtained at specific 
points in the blast wave (usually those that offer the best 
signal to noise ratio in Ha), and may not be represen- 
tative of the dynamics of the entire SNR. Furthermore, 
iHeng fc McCra v (2007) recently expressed concern that 
the current models might be underestimating the shock 
velocities from the Ha observations. For these reasons, 
the values of ufs listed in Table [3] need to be considered 
with caution. 

We have used the average ionization timescale (n e t) of 
Fe and Si to represent the ionization state of these el- 
ements in the shocked ejecta of our SNRs. In order to 
assemble a consistent set of values and constrain the un- 
certainties in the best possible way, we have determined 
the ionization timcscales listed in Table 2] from archival 
Chandra observations of our SNR sample. The methods 
and techniques that we have used for this purpose are 
detailed in Appendix [XJ 

5.2. Forward Shock Dynamics: Models vs. Observations 

Before attempting a comparison between the dynamics 
of the FS in our SNR models and the observed values, we 
must emphasize that our calculations do not include the 
effect of efficient shock acceleration of cosmic rays. This 
process is known t o take part in many young SNRs, and 
recent theoretical (lEllison et al.ll2004f) and o bservational 
pecourchelle et al.ll2000l: I Warren et al.ll2 005) works have 
shown that it can have a noticeable impact on the inter- 
nal structure and X-ray emission of the SNR, specially 
in the shocked ISM. The impact on the dynamics of the 
FS, however, is modest: at an age of 5000 yr, the values 
of Rfs and ufs m a SNR that is accelerating cosmic 
rays efficiently are only ~ 15% lower than in a test par- 
ticle case (see Fig. 3 in lEllison et~aTll2004[ ). This kind of 
deviation is inconsequential for the order of magnitude 
comparisons that we carry out in the present Section. 

The values of Rfs and Ufs m our SNR models are 
plotted alongside the observations in Figures [5] and [71 In 
Figure [51 we display the values obtained for the interac- 
tion between the explosion model PDDe and the CSM 
profiles HP3 and L2 from Figure [31 together with a set 
of SNR models generated with a uniform ISM. For the 
uniform ISM interaction, we have represented the param- 
eter space obtained by varying pism as a dashed region 
between two limiting cases, pism — 5 x 10~ 25 gcm~ 3 
(top plots) and pism = 5 x 10~ 24 gcm~ 3 (bottom plots), 
which encompass the typical conditions in the warm 
atomic phase of the ISM. The SNR models evolving in- 
side the HP3 and L2 CSM profiles are clearly unable to 
reproduce the observed values of Rfs, but most of the 
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data points cluster nicely around the uniform ISM mod- 
els. The only exception is SN 1006, whose large size is 
compatible with both a very low density ISM and the 
CSM cavity from model HP3. The situation is some- 
what different for the FS velocities. While it is clear 
that the SNR models evolving inside the CSM cavities 
can be confidently discarded, it seems that the models 
interacting with a uniform ISM also overpredict the val- 
ues of u_fs, again with the only exception of SN 1006. 
The distribution of the data points for Kepler, Tycho, 
0509-67.5 and 0519-69.0 suggests a systematic shift of a 
factor ~ 2 between the values of u fs predicted by mod- 
eling the Ha emission of the FS and the models with an 
ISM interaction. This could be th e hallmark of the ef- 
fect discussed bv lHeng fc McCravl ((2007) , or a selection 
effect towards localizations in the FS with higher densi- 
ties imposed by the statistics of the Ha observations, but 
more complex explanations cannot be discarded without 
a detailed analysis. 

In Figure [7J we repeat the comparison, this time for 
SNRs evolving into the cavities obtained with the mod- 
ified values of u f- From the point of view of the FS 
dynamics alone, the CSM structures produced by out- 
flows in the slow regime are hard to distinguish from a 
uniform ISM interaction. This is specially true for the 
collapsed cavities of models L2ule6 and L2ule7 and the 
small RB of model HP3ule6. The energy-driven cavities 
of models HP3u2e7 and L2u2e7, on the other hand, are 
hard to reconcile with the observations of all SNRs except 
SN 1006. The size of these two cavities is coincidentally 
similar to the radius of SN 1006, so that the FS is hitting 
the shell at R c for a SNR age of ~ 1000 yr. This results 
in a good match to the present value of ufs in SN 1006, 
but in a regime where the FS is undergoing a very rapid 
deceleration as it overcomes the radiatively cooled shell 
(dash-dotted plots in Figure UJ. If this was indeed the 
situation, one would expect to see a radiative FS with a 
very low expansion parame ter in SN 1006, inste ad of the 
observed nonradiative FS (I Winkler et al.ll2003l) wit h an 
expansion parameter of ~ 0.48 ( Moffett et al.lll993f ). 

We conclude this section noting that there is evi- 
dence for some kind of CSM interactio n in at least two 
of the objects in our sample, Kepler (Bandicra 119871 : 

"20031) . 



Borkowski et alJ 1199211 and N103B (jLewis et al.l 



The fact that this is not clearly revealed by the values 
of Rfs and ufs listed in Table [3] emphasizes the impor- 
tance that hydrodynamic models of individual SNRs will 
have in order to take these comparisons to the next level 
of detail. 

5.3. X-ray Emission from the Shocked Ejecta: Models 
vs. Observations 

In Figures [5] and [H we compare the values of (n e t) 
for Fe and Si listed in Table [4] to the emission measure 
av eraged ionization tim escales ((r) as defined in § 4.2 
of iBadenes et al.l 12003) in our models. These compar- 
isons have to be of a more qualitative sort than those 
based on the FS dynamics, both due to issues that affect 
the spectral models that we use to fit the X-ray obser- 
vations (see Appendix [AJ and to the limitations of our 
one-dimensio nal hydrodynamic ca lculations (see discus- 
sion in § 8 of IBadenes et al.l 120061 ). In spite of this, we 
find that qualitative comparisons of ionization timescales 
are highly informative, and more than sufficient for our 



present needs. 

In Figure [51 we present the temporal evolution of (r) 
for Si and Fe in the SNR models obtained from the in- 
teraction of model PDDe with the CSM profiles HP3 
and L2, and with a uniform ISM. For the uniform ISM 
case, we have represented the variation of pism between 
Pism — 5xl0~ 25 gcm -3 and pis m = 5xl0~ 24 gcm~ 3 as 
in Figures [6] and We have scaled (r) w ith pism using 
the ap proximate relations given in § 2.3 of Badcncs et al. 
(2005). Once again, the data points cluster around the 
models with a uniform ISM interaction, albeit with more 
dispersion than in the Rfs plots. The inconsistency be- 
tween the SNR models evolving inside low-density cavi- 
ties and the observations is even more dramatic than for 
the FS dynamics. There is simply not enough mass in the 
CSM models HP3 and L2 that the SN ejecta can react 
to in order to reach the ionization timescales observed in 
Type la SNRs. 

In Figure [51 we repeat the comparison using the CSM 
profiles obtained with the modified values of u f . As ex- 
pected, the CSM profiles generated by outflows in the 
slow regime lead to a stronger interaction with the SN 
ejecta and higher ionization timescales for Fe and Si. 
This is specially true for models HP3ule6 and L2ule6, 
which have prominent shells of radiatively cooled mate- 
rial close to the SN progenitor. In these models, some 
ejecta layers undergo catastrophic cooling early on the 
evolution of the SNR. Models HP3ule7 and L2ule7 do 
not lead to such extreme interaction, and are in gen- 
eral hard to distinguish from a uniform ISM. As soon as 
the outflows enter the fast regime, however, the situa- 
tion changes radically. Only the shocked ejecta emission 
of SN 1006, with its low Si n e t and absence of Fe, is 
compatible with the CSM cavities HP3u2e7 and L2u2e7, 
but these models have already been discarded in § 15. 21 for 
other reasons. 

6. DISCUSSION 

In § [3j we have seen that the accretion wind outflows 
invoked by current Type la progenitor models in the SD 
channel will excavate large cavities in the ISM. In § [SJ we 
have shown that the fundamental properties of the seven 
Type la SNRs in our sample are incompatible with SNR 
models that expand into such large cavities. Can these 
two results be connected with enough confidence to draw 
conclusions on Type la SN progenitor models based on 
the SNR observations? In this Section, we review the 
potential issues that might affect the initial models, the 
simulations, and the observational sample, and we briefly 
discuss the implications of our results for the observa- 
tions of Type la SNe. 

6.1. Initial Models: The True Nature of Accretion 
Winds. 

The outflow models presented in Figure [1] and Table [1] 
make two important approximations: they are isotropic 
(i.e., spherically symmetric) in space and continuous in 
time. In a more realistic scenario, the presence of the 
donor star and the accretion dis k, and the spin-up of 
the W D to high angular velocities (|Yoon fc Langerll2004l 
|2005t ) should introduce some sort of bipolarity in the out- 
flows. As we have seen in § 12.2.21 the observations of 
Type la SN progenitor candidates show some evidence 
of this bipolarity. Continuous bipolar outflows may lead 
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to a bip olar CSM structure s imilar to some planetary 
nebulae (jBalick fc Frank! |2002| ). but this does not seem 
likely for the objects in our sample. Only moderate one- 
sided asymmetries are apparent in Kepler and N103B, 
and most other objects (in particular, Tycho, SN 1006, 
and 0509-67.5) are nearly spherical. 

Setting aside the spatial structure of the outflows from 
Type la SN progenitors, their temporal evolution will 
not be as simple as the plots shown in Figure [TJ In or- 
der to reach the values of M that allow for the growth 
of the WD, the mass transfer rate has to cross the in- 
stability region below M sta Ue, wh ere more or less fre - 
quent nova outbursts are expected. lLanger et al.l (2000), 
for instance, found a long (~ 10 6 yr) switch-on phase of 
the mass transfer in all their binary systems that was 
dominated by instabilities. Nova-like outbursts could 
also appear during the mass-conservative phase between 
the c essation of the accret ion wind and the SN explo- 
sion (jHachisu et al.l [199 9^. Judging from the obser- 
vations of RXJ0513. 9-6951 and V Sagittae discussed in 
§ 12.2.21 a steady accretion wind might be impossible or 
happen only in a few cases. If the episodic outflows in 
Type la progenitors have periodicities of ~ 100 days, 
as in these two systems, several million cycles should 
take place through the accretion phase, and the outflows 
might smear out to become indistinguishable from a con- 
tinuous wind. In the case of recurrent novae, the indi- 
vidual outbursts will be more spaced out, leading to a 
complex CSM structure, with several shells that brighten 
up as they are overtaken by the FS. Again, there is no 
evidence for this in any of the objects that we have stud- 
ied, but more detailed hydrodynamic calculations would 
be necessary to confirm or discard any given scenario. 

6.2. Simulations: From Outflows to Cavities. 

There are two ways to prevent Type la progenitor out- 
flows from leaving large cavities in the CSM: either the 
value of u f is decreased below the critical limit u cr or the 
cavities are somehow destroyed before the SN explodes. 
As we have shown in § [5j outflows in the slow regime 
(u Q f < 100 km s" 1 ) lead to CSM structures that are hard 
to distinguish from a uniform ISM interaction. However, 
outflows emanating from the WD surface should have 
much higher velocities, of the order of the escape velocity 
(~ 1000 km s -1 , a value that is supported by the obser- 
vations cited in § I2.2.2[) . Some of the kinetic energy of 
the outflows will be lost in overcoming the potential well 
of the system, but it is unlikely that this can reduce u a f 
by an order of magnitude. Furthermore, the presence of 
slow outflows around the progenitor systems is difficult 
to reconcile with the lack of CSM detections in Type la 
SNe (see §[13). 

If large cavities do form, it is hard to find a non- 
catastrophic process that could destroy them before tsN- 
Thermal conduction, for instance, which must be ac- 
tive to some extent at the interface between the hot 
shocked outflow and the radiatively cooled shell, could 
not wipe out cavities with a radius of ~ 10 20 cm, even 
over ti mescales of ~ 10 6 yr. In the example presented in 
§ IV of IWeaver et al.1 (|1977l ). which is based on an inter- 
stellar bubble very similar to our models, the conduction 
front behind the cool shell only spans a small fraction of 
the cavity radius. The density inside the cavity increases 



due to evaporation of material from the shell, but it re- 
mains orders of magnitude below the surrounding ISM. 

A final possibility is that the binary system forms a 
cavity, and then leaves it behind as it moves with respect 
to the ISM. Large speeds are not expected in SD Type 
la progenitors, but even at the mean random velocity of 
~ 20kms~ 1 (character istic of the solar neighborhood, 
iDehnen fc Binne"vlll998| ) a star can move ~ 10 20 cm in 
2 x 10 6 yr. Still, only moving systems with an extended 
mass conservative phase before the SN explosion could 
leave their cavities completely behind. In moving pro- 
genitors with active outflows at tsjv, the explosion would 
happen off-center, but the cavity would still be there, and 
therefore the impact on the FS dynamics and the X-ray 
emission of the SNR would still be strong (| Weaver et al.l 
119771: iRozvczka et al.lfl993h . 

6.3. Observations: Other Type la SNRs? 

We have omitted from our study several SNRs whose 
X-ray spectra suggest a Type la origin, but whose ages 
can only be roughly estimated from the FS dynam- 
ics. Exam ples include G337.2-0 7 and G299.2-2 .9 in 
the Galaxy (|Rakowski et alJl2006t|Park et al.ll2006l ). and 
DEM L71, DEM L238, and DEM L249 in the LMC 
([Hughes et al.1 120031 : iBorkowski et all 120061 ) . These are 
rather old objects, and nothing in their morphology or 
X-ray emission indicates that they might be expanding 
inside a large wind-blown cavity, but clues for other kinds 
of CSM interaction might be revealed by detailed models. 

It is also possible that young Type la SNRs expanding 
inside large cavities do exist, but have not been identified 
yet. These objects would have a large radius and a faint 
X-ray emission, so they would be hard to classify as Type 
la SNRs. One such object might be the Galactic SNR 
RCW 86, whic h has a large size (<xfs = 42') a nd Fe- 
rich SN ejecta ([Vink et all 1 19971 : IRho et alJl2002f ). This 
SNR has been ass ociated with the supernova of AD 185 
(|Vink et al.ll2006l ). but this association remains contro- 
versial, and its Type la origin is uncertain. If we ignore 
th ese concerns, and tak e the value of D = 2.8 kpc given 
bv iRosado et alj (|1996f) . we find R FS = 1.1 x 10 20 cm 
at an age of 1821 yr, about 30 % larger than model 
PDDe+L2, but certainly within the correct order of mag- 
nitude. Recent va lues for the FS velocity (~ 2700 km s , 
IVink et ail l2006h and the i onization timesc ale of the 
shocked ejecta (~ 10 9 cm" 3 s, IRho et al.H"2002f ) are much 
closer to our cavity models PDDe+HP3 and PDDc+L2 
than t hose of the SNRs in our sample. In fact. I Vink et al.l 
(|1997h already suggested that RCW 86 might be evolving 
inside a low-density cavity based on early ASCA observa- 
tions. If the explosion type and the association with SN 
185 are confirmed for RCW 86, it would be very interest- 
ing to explore the possibility of an accretion wind outflow 
as the origin of the CSM structure. Another Galactic 
SNR with Fe-rich ejec ta that might be expanding into a 
CSM cavity is W49B ([Hwang et alJl2000t iKeohane et all 
l2007h . but in this case there is no reliable age estimate 
that can be used to constrain the dynamics. 

6.4. Implications of accretion wind outflows for the 
observations of Type la SNe 

Although this paper focuses on the impact that the 
outflows from Type la progenitors have on the dynamics 
and X-ray emission of SNRs, here we will briefly discuss 
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the implications for observations of Type la SNe. The 
lack of prompt emission at X-ray and ra dio wavelengths 
(|Panagia et al.|[2~006t llmmler et al.ll2006l and other refer- 
ences listed in §[!]), as well as the absence of low velocity 
(narrow) li nes from H in the ea rly or late spectra of nor- 
mal events (|Mattila et al.ll2005l ) suggest that the amount 
of circumstellar material in the immediate vicin ity of 
Type la progenitors must be low (|Eck et al.l 12002). The 
upper limits on the radio and X-ray fluxes have been used 
to derive estimates for the mass loss rate of the progeni- 
tor system outflows M /, but these observations can only 

constrain the quotient M o f/u f- If the outflows are fast, 
this technique does not provide much information about 
M f- the radio fluxes of iPanagia et a l. (2006) translate 
into M of < 3 x lO^Mgyr" 1 for u of = 1000 km s"\ 
which is compatible with all but the most extreme of ac- 
cretion winds. For slower outflows, the upper limits on 
M f go down li nearly with u a f. An interesting possibil- 
ity, suggested bv lWood-Vasev fe Sokoloskl (|2006l ) , is that 
episodic outflows just before the explosion might clear 
out a more or less extended region around the progenitor 
system, thus explaining the lack of prompt emission and 
H-rich material. However, we have shown that this re- 
gion cannot be so large that the dynamic evolution of the 
SNR is severely affected, at least in the seven objects that 
we have studied. At the same time, evidence for large 
CSM structures around some Type la progenitors has ap- 
peared in the form of the light echoes detected in three 
objects: SN 1991T, SN1995E a nd SN 1998bu. Th ese 
echoes have been interpreted by iQuinn et alJ (|2006l ) as 
reflected light from detached shells or sheets of dust at 
distances ranging between 50 and 20 pc from th e SN. 
If this interpretation is correct (see iPatatJ 120051 . for a 
complete discussion of the ambiguities involved) , the de- 
tached structures in these three objects might be the ra- 
diatively cooled shells around accretion wind-blown bub- 
bles. 

7. CONCLUSIONS 

In this paper, we have explored the relationship be- 
tween the SD progenitor systems of Type la SNe and 
their surroundings, focusing on the effects that the out- 
flows from these systems would have on the structure of 
the CSM and on the dynamics and X-ray emission of 
the SNRs that evolve after the explosion. We have seen 
that the optically thick accretion winds from the WD 
surface invoked by current Type la progenitor models 
will excavate large cavities in the ISM, provided that the 
outflows stay in the fast regime (u a f > 200 km s -1 ). We 
have shown that the fundamental properties of the seven 
young Type la SNRs in our sample (SN 1885, Kepler, 
Tycho, SN 1006, 0509-67.5, 0519-69.0, and N103B) are 
incompatible with SNR models that expand into such 
large cavities. In fact, we found that all these objects 
can be explained by a uniform ISM interaction at the 
level of detail allowed by our simulations. At the same 
time, we cannot discard the existence of a population of 
young Type la SNRs expanding into low-density cavities 
excavated by accretion wind outflows. The properties 



of these objects (large size and low brightness) would 
make them hard to detect and, if detected, very diffi- 
cult to identify as the remnants of Type la explosions. 
The Galactic SNR RCW 86 might be such an object, but 
more detailed work is needed to confirm this intriguing 
possibility. 

The presence of outflows from the progenitors of the 
seven SNRs that we have examined cannot be completely 
discarded, but the properties of these outflows (mass-loss 
rates, timescales and velocities) must be very different 
from those of the accretion winds described by HKN96. 
It remains to be seen whether outflows with substantially 
lower values of M a f or u f would be capable of limit- 
ing the effective accretion rate of the WD and prevent 
a common envelope phase in single degenerate Type la 
progenitors. Another possibility is that accretion wind 
outflows are not always present in these systems - for 
instance, if the rotation of the WD increases the accre- 
tion efficiency to the point that only a small fraction 
of systems pres e nt sig nificant outflows, as proposed by 
lYoon fe Langen Ip005fi . 

We conclude with a reminder that the viability of SD 
systems as Type la progenitors has not been proved yet. 
Several studies have pointed out the difficulties in driv- 
ing a H- accreting WD clo s e enough to the Chandrasek har 
mass (ICassisi et ; al.lll998t iPiersanti et al.lfl999L l2000D. If 
the star found bv lRuiz-Lapuente et al.l ( 2004 ) in the Ty- 
cho SNR is indeed the runaway companion of the WD 
that exploded in 1572, this means that at least one of 
the objects in our sample must have had a SD progen- 
itor. The obvious alternative (DD systems) is plagued 
by its own problems, and it is still uncle ar whether WD 
mergers can lead to health y explosions (|Segretain et al.l 
Il997t iGuerrero et al.ll2004h . The fact remains that we 
lack a clear picture of the binary evolution leading to 
Type la SNe. Studies of the CSM structure in young 
Type la SNRs should provide valuable insights into this 
problem. 
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Fig. 1. — Temporal evolution of the mass loss rates in the outflow models for Type la SN progenitor systems detailed in Tablefl] Left: 
models HI, LV1, LI, and L2. Right: models HP1, HP2, and HP3. The symbols mark the time of the SN explosion, tsN- 

Yoon, S.-C, Langer, N., & van der Sluys, M. 2004b, A&A, 425, 
207 



APPENDIX 

IONIZATION TIMESCALES FOR THE SHOCKED EJECTA IN TYPE IA SUPERNOVA REMNANTS 

Together with the electron temperature T e , the ionization timescale or fluence (n e t) is a fundamental quantity used 
to characterize the thermal emission in nonequilibrium ionization (NEI) plasmas. These two parameters are the most 
common end product of fitting the thermal component of the X-ray emission in SNRs. However, the values of T e 
and n e t obtained in a particular spectral fi t are model dependent , and can be affected b y uncertainties in the atomic 
data and other factors (for discussions, see iBorkowski et ail 1200 lURakowski et alJl2~006f ). The most popular tools for 
spectral analysis of the thermal emission in SNRs are the plane-parallel shock models with adjustable abundances, but 
even these sophisticated spectral models contain fundamental simplifications whose impact on the fitted parameters is 
unclear, like the assumption of a homogeneous chemical composition or the plane shock geometry In general, it can 
be said that the systematic uncertainties associated with plane-parallel shock fits are much larger (and much harder 
to estimate) than the statistical ones. 

In order to minimize these sources of error and determine values of n e t in different SNRs that can be compared to each 
other (and to the models) with some confidence, we have downloaded archival Chandra observations for the objects 
in our sample, and we have fitted them with the same procedure. We used the plane-parallel shock models described 
in iHughes et a l. (2000) for the thermal component and a power law to describe the nonthermal continuum, both of 
which were absorbed by an intervening column of Galactic interstellar matter assumed to be of solar composition. 
During the fits, the Galactic column density was fixed to a value determined from the entire Chandra spectrum. To 
minimize the errors related to atomic data, we concentrated on the K-shell lines from the most abundant elements 
(spectra above 1.5 keV). This has the advantage of reducing the impact of uncertainties in the absorption as well. In 
fits to the 1.5-5 keV band ('Si component') we included emission from Si, S, Ar, and Ca, all with the same T e and 
n e t. The relative abundances of these species were allowed to vary during the fit, as were the index and normaliztaion 
of the power-law. Independent fits were done to the 5-10 keV band around the Fe K blend ('Fe component'). Only Fe 
and the power-law continuum (plus fixed absorption) were included in the Fe fits. The best-fit models for each SNR 
in both bands can be seen in Figure 1101 

The relationship between the two fundamental quantities n e t and T e in these fits is of particular interest. Plane- 
parallel shock models implicitly rely on the thermal continuum to constrain T e and the line emission to constrain n e t. 
Unfortunately, the X-ray continuum is often dominated by nonthermal emission (the power-law in our fits), and the 
value of T e is largely unconstrained. We have taken this situation into account by producing fits at fixed values of T e : 1 
and 10 keV for the Si component, 5 and 10 keV for the Fe component. These values represent reasonable maxima and 
minima for shock-heated plasma in young SNRs, and the best-fit value for T e (models plotted in Figure [TO]) is bounded 
by them in all cases. The range of ionization timescales listed in Table |4] corresponds to these extrema. We have also 
subtracted 0.3 dex from the fitted values to account for the fact that the plane-parallel shock models contemplate a 
distribution of ionization timescales between and n e t, so the average ((n e t)) is one half of the fitted value. 

A rough idea of the level of agreement between our hydrodynamic models and the spectral fits to the observations 
can be obtained by comparing the fitted values of log((n e t)) in the Tycho SNR from Table [4] (10.22-10.99 for Si , 
9.72-9.78 for Fe) with the values of log({r)) in the best model for the ejecta emission found bv lBadenes etafl (p00l ) 
(10.12 for Si, 9.69 for Fe). These small deviations support the validity of the comparisons made in § 15.31 

We conclude with an unrelated (but important) remark on the fitted ionization timescales listed in Table HI In all 
the Type la SNRs we have examined, the (n e t) of Si is significantly higher than that of Fe, sometimes by more than 
an order of magnitude. This constitutes strong evidence that the ejecta of Type la SNe must be stratified to some 
degree (with most of the Fe interior to most o f the Si), favoring de layed detonation models over deflagration models 
with well-mixed ejecta (see discussion in § 3 of lBadenes et^[2005h . 
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Fig. 2. — Critical outflow velocity u cr as a function of outflow velocity u f. The striped region represents the values of u cr for 
typical accretion wind outflows expanding into the warm, atomic phase of the ISM. The limiting curves are for M f = 10 — 6 Mq yr —1 , 
PISM = 5 X 10~ 24 gcm -3 (top curve) and M a f = 10 -7 Mq yr -1 , Pism = 5 X 10~ 25 gcm~ 3 (bottom curve). The dotted line has a slope 
of 1. 
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generated by outflows with modified velocities. Top panels: models HP3ule6 (dotted line), HP3ule7 (dashed line), HP3u2e7 (dash-dotted 
line). Bottom panels: models L2ule6 (dotted line), L2ule7 (dashed line), L2u2e7 (dash-dotted line). 
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Fig. 8. — Temporal evolution of the emission measure averaged ionization timescale (t) for Si (left panel) and Fe (right panel) in the 
shocked ejecta of model PDDe interacting with a constant density ISM (solid lines), and with the CSM profiles HP3 (dotted line) and L2 
(dashed line). The two solid plots correspond to an interaction with Pism = 5 X 10 — 25 gem -3 (top plots) and Pism = 5 X 10 -24 gem -3 
(bottom plots). The striped region between the solid lines spans the parameter space of pism between these two values. The data from 
Table [4] are overlaid on the plots. 
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Fig. 9. — The same as Figure[8] but for the explosion model PDDe interacting with a uniform ISM (solid lines), and the CSM profiles 
generated by outflows with modified velocities. Top panels: models HP3ule6 (dotted line), HP3ule7 (dashed line), HP3u2e7 (dash-dotted 
line). Bottom panels: models L2ule6 (dotted line), L2ule7 (dashed line), L2u2e7 (dash-dotted line). 




Fig. 10. — Plane-parallel shock fits to the Chandra spectra of the SNRs in the sample. The K-shell emission blends from Si, S, Ar, and 
Ca are on the left panel, and the Fe K blend is on the right panel. The spectra are roughly ordered from highest (top) to lowest (bottom) 
ionization timescale. 
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TABLE 1 

Outflow Models for Type Ia Progenitor Systems 



Model M f tsN Binary System Parameters 

Name (M ) (yr) M WDfi (M G ) M D>0 (M ) P (days) 



Reference 



HI 
LV1 
HP1 
HP2 
HP3 
LI 
L2 



0.15 
0.50 
0.24 
0.80 
0.50 
0.40 
0.64 



5.0 x 10 5 
1.8 x 10 6 
2.0 x 10 6 
1.4 x 10 6 
6.0 x 10 5 
2.0 x 10 6 
2.0 x 10 6 



1.0 
1.0 
0.75 
0.8 
1.0 
1.0 
0.8 



2.0 


2.0 


1 (Fig. 7) 


2.5 


1.6 


2 (Fig. 1) 


2.0 


1.58 


3 (Fig. la) 


2.2 


2.50 


3 (Fig. lc) 


2.4 


3.98 


3 (Fig. lc) 


2.3 


1.74 


4 (Model 2, Fig.7) 


2.1 


1.53 


4,5 (Model 31, Fig. 36 in ref. 5) 



R eferences. — (l): | Hachisu et a ll lll999b l) 
(4): ILanger et~aT1 l l200oh ; (5): IDeutschmannl 11199^ 



Li & van den Houvcl (1997); (3): Han & Podsiadlowski (2004); 



Note. — Outflow model parameters: M f. total mass lost by the system; tsN' time between the onset of the 
final mass transfer episode and the SN explosion. Binary system parameters: Mwd,o'- initial mass of the WD; 
q: initial mass of the donor; Pq: initial period. 



TABLE 2 

Structure of the cavities produced by Type Ia Progenitor Outflows 



Outflow Parameters Cavity Parameters 







PlSM 




Outflow 


Ruof 




Rc 




Stage 


Model 


(cms 1 ) 


(gcm~ 3 ) 


(K) 


Regime 


(cm) 


(M Q ) 


(cm) 


(M Q ) 


at tsN 


HI 


10 s 


lfJ -24 


10 4 


F 






5.25 x 10 19 


3.05 x 10 2 


APCB 


LV1 


10 s 


1Q -24 


If) 4 


F 






9.75 x 10 19 


1.95 x 10 3 


APCB 


HP1 


10 8 


1Q -24 


If) 4 


F 






8.30 x 10 19 


1.20 x 10 3 


APCB 


HP2 


10 8 


1Q -24 


10 4 


F 






1.06 x 10 20 


2.51 x 10 3 


APCB 


HP3 


10 s 


10 -24 


10 4 


F 


1.40 X 10 19 


8.00 x 10~ 3 


6.68 x 10 19 


7.68 x 10 2 


ABROS 


LI 


10 s 


1Q -24 


10 4 


F 


1.17 x 10 19 


1.45 x 10~ 3 


9.23 x 10 19 


1.66 x 10 3 


APCB 


L2 


10 s 


1Q -24 


10 4 


F 






1.10 x 10 20 


2.80 x 10 3 


APCB 


HP31owp 


10* 


5 x lO" 25 


5 x 10 3 


F 


1.80 x 10 19 


1.03 x 10~ 3 


7.82 x 10 19 


5.97 x 10 2 


ABROS 


L21owp 


10 8 


5 x 10~ 25 


5 x 10 3 


F 






1.35 x 10 20 


2.59 x 10 3 


APCB 


HP3ule6 


10° 


1Q -24 


10 4 


S 


3.75 x 10 18 


0.28 


4.89 x 10 18 


1.64 


RB 


HP3ule7 


10 7 


1Q -24 


If) 4 


S 


1.18 x 10 19 


8.10 x 10~ 2 


1.70 x 10 19 


0.65 


PRB 


HP3u2e7 


2 x 10 7 


1Q -24 


If) 4 


F 


1.14 x 10 19 


3.75 x 10~ 2 


3.13 x 10 19 


1.44 x 10 2 


ABROS 


L2ulc6 


10 6 


10 -24 


10 4 


s 






9.80 x 10 17 a 


2.45 a 


RB (C) 


L2ulc7 


10 7 


1Q -24 


10 4 


s 










PRB (C) 


L2u2c7 


2 x 10 7 


1Q -24 


10 4 


F 






3.13 x 10 19 


52.5 


APCB 


Note. — 


Outflow re| 


$imc: fast (F] 


or slow (S). Cavity parameters: R r 


of. radius of the unshocked outflow; M uo f. 


mass of the 



unshocked outflow; R c : radius of the CD (for models in the fast outflow regime, this is equal to the inner radius of the low-density 
cavity); M s h' mass in the radiatively cooled shell (for cavities that are not pressure confined, M s h is the mass c o ntained between 
R c and the outer supersonic blast wave). For the cavity stage at tsN, we use the terminology of lKoo & M cKcc (1992b). APCB: 
adiabatic pressure-confined bubble; ABROS: adiabatic bubble with a radiative outer shock; RB: radiative bubble; PRB: partially 
radiative bubble. Cavities marked with '(C)' have collapsed at tsN- 
a Parameters for the relic shell of the collapsed bubble. 



TABLE 3 

Type Ia Supernova Remnants - Forward Shock Dynamics 



Remnant 


Age 


D 


«FS 


Rfs 


ups 


Name 


(yr) 


(kpc) 


(arcmin) 


(cm) 


(cms -1 ) 



SN1885 


121 a 


785 ±30 (1) b 


6.7 x 10~ 3 (1) 


(4.7 ±0.2) x 10 18 




Kepler 


402 a 


4.8 ± 1.4 (2) 


1.7 


(7.3 ±2.1) x 10 18 


2.0 x 10 8 - 2.5 x 10 8 (3) 


Tycho 


434 a 


1.5 - 3.1 (4) 


4.0 


(8.3 ± 2.9) x 10 18 


1.9 x 10 s - 2.3 x 10 8 (5) 


SN1006 


1000 a 


1.4-2.8 (4) c 


15.0 


(2.8 ± 0.9) x 10 19 


2.8 x 10 s - 3.0 x 10 s (3) 


0509-67.5 


400 ± 120 (7) d 


50 ± 1 (8) c 


0.24 f 


1.1 x 10 19 


> 2.0 x 10 8 (5) 


0519-69.0 


600 ± 200 (7) d 


50 ± 1 (8) c 


0.26 f 


1.2 x 10 19 


1.0 x 10 8 - 1.9 x 10 8 (5) 


N103B 


860 ±400 (7) d ' g 


50 ±1 (8) e 


0.23 f 


1.0 x 10 19 





References. — (1) | Fesen et"aT1 j2U06h: (2) iRevnoso fe Goss) (fl999h : (3) ISollerman et al.l d2003tt : f4) ISmith et alj (fl99lT) : (5) 
Ghavamian et al 

] (I200lh : (7)\ Rest et al 1 (120051 ): f8l lAlvesl ll20ol i 
Note. — SNR parameters: D: distance to the SNR; qfs: angular radius of the FS; Rfs'- FS radius (Rfs — 3.09 x 
10 21 L>( QFS /60)(7r/180)); u F s- FS velocity. 

a Historical SNR. b Distance to M31. c For the distance to SN1006, we use the more conservative estimate of lSmith et al.| d!991fh 
which includes systematic as well as statistical uncertainties, instead of the newer results of Winkler et al. (2003).^ Age estimate 
from light echoes. Distance to the LMC/ For the LMC SNRs, the radii have been determined by the outermost extent of the 
X-ray emission in the 0.5-4 kcV band, sub tracting 0.5" due to the smearing of the rim by the Chandra PSF. g The uncertainty in 
the age of N103B arbitrarily set to 400 yr. iRest et al ctermincd an age of 860 yr, but could not constrain the uncertainty 

in the measurement. 
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TABLE 4 
Type Ia Supernova Remnants - 
Ionization Timescales in the Shocked 
Ejecta 



Remnant 


l°f 


5 {net) ( 


cm 3 s) 


Name 


Si 


Fe 


Kepler 


10.08 - 


10.24 


9.85 - 9.92 


Tycho 


10.23 - 


10.99 


9.72 - 9.78 


SN1006 


9.49 - 


9.60 


a 


0509-67.5 


9.80 - 


9.82 


< 9.7 


0519-69.0 


10.50 - 


11.62 


9.90 - 9.95 


N103B 


10.64 - 


11.94 


10.62 - 10.69 



a SN 1006 has no Fc K emission in its Chandra 
X-ray spectrum (Hughes et al., in preparation). 



